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We report on a general theory for analyzing quantum transport through devices in the Metal- 
QD-Metal configuration where QD is a quantum dot or the device scattering region which contains 
Rashba spin-orbital and electron-electron interactions. The metal leads may or may not be fer- 
romagnetic, they are assumed to weakly couple to the QD region. Our theory is formulated by 
second quantizing the Rashba spin-orbital interaction in spectral space (instead of real space), and 
quantum transport is then analyzed within the Keldysh nonequilibrium Green's function formalism. 
The Rashba interaction causes two main effects to the Hamiltonian: (i) it gives rise to an extra 
spin-dependent phase factor in the coupling matrix elements between the leads and the QD; (ii) 
it gives rise to an inter-level spin-flip term but forbids any intra-level spin-flips. Our formalism 
provides a starting point for analyzing many quantum transport issues where spin-orbital effects are 
important. As an example, we investigate transport properties of a Aharnov-Bohm ring in which a 
QD having Rashba spin-orbital and e-e interactions is located in one arm of the ring. A substantial 
spin-polarized conductance or current emerges in this device due to a combined effect of a magnetic 
flux and the Rashba interaction. The direction and strength of the spin-polarization are shown to 
be controllable by both the magnetic flux and a gate voltage. 

PACS numbers: 72.25.Dc, 73.23.-b, 85.75.-d, 73.40.Sx 
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I. INTRODUCTION 

The spin-orbit (SO) interaction in semiconductors has 
attracted great attention in recent yearsii2i2*^iS*2iL2i2ii2iii 

as it plays a very interesting role for the emerging field 
of semiconductor spintronics. SO interaction can couple 
the spin degree of freedom of an electron to its orbital 
motion and vice versa thereby giving a useful handle for 
manipulating and controlling the electron spin by exter- 
nal electric fields or gate voltages. SO is an intrinsic 
interaction having its origin from relativistic effects, but 
it is believed to be substantial in some semiconductors. 
More than ten years ago, Datta and Das theoretically 
analyzed the possibility of a spin-transistor which works 
due to Rashba SO interaction that induces spin preces- 
sions in a semiconductor— with ferromagnetic leads. More 
recently, Murakami et.al. theoretically predicted that a 
substantial amount of dissipationless quantum spin cur- 
rent could be generated by a co-action of electric field 
and SO interaction^^. Shen et.al. found a resonant spin 
Hall conductance in a two-dimensional (2D) system with 
Rashba SO interaction under a perpendicular magnetic 
field^. There are also many other works on related issues 
where SO interaction plays a central roleiSiiiiiiiSiiSiil, 
and this research direction is expanding in a very rapid 
pace due to its possible application to spintronics. 

A semiconductor spintronic device is likely to be 
based on mesoscopic and nanostructures where electron- 
electron (e-e) interactions may be strong. Hence it is 



desirable to formulate a general quantum transport the- 
ory which can handle SO, e-e, and other interactions 
for systems in the Metal-QD-Metal configuration. Here 
the "QD" indicates a quantum dot or the device scat- 
tering region where the various interactions exist, while 
the "Metals" are device leads which extend to electron 
reservoirs far away. The metal leads may or may not be 
ferromagnetic, but are weakly coupled to the QD region. 
In almost all previous theoretical work, the SO interac- 
tions are represented by a real space Hamiltonian where 
e-e interactions and strong correlations are usually ne- 
glected. Indeed, it is rather difficult, if not impossible, to 
handle SO together with e-e correlations and other inter- 
actions in real space for transport problems. In contrast, 
the most powerful and general theoretical technique for 
quantum transport in meso- and nano-scopic systems is 
the Keldysh nonequilibrium Green's function (NEGF) 
formalismi^. NEGF can handle many-body correlations 
and interactions in an unified fashion and is a well es- 
tablished formalism-'^^. NEGF is typically formulated in 
momentum space or other spectral space for theoretical 
and numerical analysis. This means that all interactions 
need to be formulated in the spectral space. In other 
words, in order to conveniently apply NEGF theory, one 
needs to write the SO interactions in a spectral space 
with second quantization. To the best of our knowledge, 
we are not aware of a derivation of such a second quan- 
tized form for SO interaction. 

It is the purpose of this paper to report a general quan- 
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turn transport theory for Metal-QD-Metal devices with 
SO and e-e interactions, based on the NEGF framework. 
We begin by presenting a pedagogical discussion of the 
general physics of SO interaction by quantizing the cor- 
responding classical action, which gives a vivid physical 
picture of SO interaction. We then second quantize the 
real space SO interaction Hamiltonian in a proper spec- 
tral space so that it can be analyzed by NEGF for quan- 
tum transport properties of Metal-QD-Metal devices. Al- 
though the derivations are general, we specialize on a 
Metal-QD-Metal device where the QD is described by 
the Anderson model plus the Rashba SO interaction, and 
the leads are ferromagnetic metal. The second quantized 
Hamiltonian can then be analyzed within NEGF and well 
established many-body theoretical methodsi^iSfliSiiSS. To 
illustrate our formalism, we investigate quantum trans- 
port properties of a Aharnov-Bohm ring where a QD 
having Rashba SO and e-e interactions sits in one arm 
of the ring. We found that a substantial spin-polarized 
conductance or current emerges in this device when a 
magnetic flux passes through the ring. In particular, its 
spin-polarized direction and strength are controllable by 
both the magnetic flux and a gate voltage, hence the 
predictions are testable experimentally. 

The paper is organized as follows. In the next sec- 
tion we discuss the many-body Hamiltonian of a general 
Metal-QD-Metal device structure, and present a peda- 
gogical discussion of the SO interactions in real space. 
We then proceed, in Section HI, to second quantize the 
Rashba SO interaction in spectral space so that the entire 
device Hamiltonian can be second quantized. This pro- 
cess is divided into several subsections and careful deriva- 
tions and discussions will be presented. A brief summary 
for the second quantized Hamiltonian in spectral space is 
given in Section II-F. In section IV, we analyze quantum 
transport properties of a modified AB ring which con- 
tains a QD in one arm of the ring, and both Rashba SO 
and e-e interactions exist in the QD. Finally, Section V 
summarizes the results of our work. 



II. HAMILTONIAN OF THE 
METAL-QD-METAL DEVICE 

In this section we discuss the Hamiltonian of a gen- 
eral device structure in the Metal-QD-Metal device con- 
figuration. By presenting a very useful pedagogical dis- 
cussion on the classical forces acting on moving charges 
and spins inside electrical and magnetic fields, we realize 
that the SO interaction originates from the force (torque) 
of electrical field on the moving spin. This allows us to 
write down the Hamiltonian for the SO interaction in real 
space, i/so(r), and in particular we derive the Rashba SO 
interaction. Of course, the derived Hsoi'r) is the same 
as that from the Dirac equation^, but the pedagogical 
discussion gives a vivid physical picture about SO inter- 
actions for quantum transport in solid state devices. In 
fact, in the literature of relativistic quantum mechanics 



such as the book of Bjorken and Drell^- , SO interaction 
has been discussed with the point of view of quantizing 
the classical force acting on the moving spin by the exter- 
nal electric field. We found such a pedagogical discussion 
in the context of solid state electronics to be very useful. 

The general device structure we consider is schemati- 
cally shown in Fig. la, where the scattering region (QD) 
is connected to the outside world by coupling to two fer- 
romagnetic metal leads (FM). The permanent magnetic 
moments of the leads are denoted by vectors where 
(3 = L,R indicate the left and right leads. The QD is 
assumed to be weakly coupled to the leads due to the po- 
tential barriers at the two Metal-QD junctions (Fig. lb). 
Inside the QD there are SO and e-e interactions, while 
these interactions are neglected in the leads. There may 
also be an external magnetic field B(r). For this device, 
the total many-body Hamiltonian can be written as 

77(ri,r2,...,rjv)=^i/.(ri)+ ^ i?/(r„r,) (1) 

where the second term is the e-e interaction Hj{ri,r2) 
and the first term is from the single particle Hamiltonian 
Hsir): 



H,{r) = -^ + y(r)+a.M(r)+a.B(r) 

2m* 

+H,o{r) . (3) 

Hg contains the usual single particle terms: the kinetic 
and potential energies, the interaction energy with the 
magnetic moment M in the ferromagnetic leads, and 
the Zeeman energy. The last term in EqQ is the SO 
interaction Hsq. Although the real space form of Hgo 
is known^^i^ from Dirac equation, in the following we 
present a pedagogical discussion about it. 

Transport in our device is about the motion of two 
entities — charge and spin, in two fields — electric and 
magnetic fields. Therefore there are a total of four actions 
due to the fields on the charge/spin: (i) the electric or 
Coulomb force on charge; (ii) the Lorentz force on mov- 
ing charge; (iii) the magnetic force on spin (Zeeman); 
and (iv) the electric force on moving spin. Of these four 
actions, (i)-(iii) are well known and familiar, but (iv) is 
much less so. Where does (iv) comes from? It comes 
due a purely relativistic effect^. Consider a spin which 
produces a magnetic field in the space surrounding it: 
if this spin is moving, by relativistic transform we ob- 
tain an electric field (in addition to the magnetic field). 
In other words, a moving spin produces an electric field. 
Conversely, if a spin is moving inside an external electric 
field, it will be subjected to an action (torque). In this 
sense, (iv) is the counter part of the Lorentz force. It has 
been shown that a moving spin with velocity v inside an 
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electric field E is subjected to a torque action with an 
interactive potential energysSSiSl 



eh 



4m 



cr • (v X E) 



(4) 



where rrif. is the electron mass, c the speed of light, and 
the electron spin is s — ^a. Of course, this is a classical 
expression. 

To quantize the classical torque energy Q), we make 
the following replacements: electric field E Vy(r)/e 
where V{y) is the potential energy of the system; the 
speed V ^jnie. where p is the momentum operator. 
The quantum mechanical correspondence of expression 
Q becomes: 
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Finally, if we consider other forms 



of the potential energy y(r), we obtain other kinds of 
SO interactions, but the essence of SO coupling is the 
interaction of the external electric field on the moving 
spins. 



III. SECOND QUANTIZATION OF THE 
DEVICE HAMILTONIAN 



In this section we second quantize the device Hamilto- 
nian (Q). The focus is to derive the second quantization 
of the Rashba SO interaction in a spectral form. 
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[ct • (p X VT^(r)) - (7 • (VV"(r) x p)] 

[ct • (p X VT^(r)) + VF(r) .(ax p)] (5) 



where a = [a^, cTy,az) is the vector of Pauli matrix. Ex- 
pression (0 is exactly the general form of the SO interac- 
tion Hamiltonian, usually derived from the Dirac equa- 
tion in the low speed limits^^. Therefore, the essence of 
the SO interaction is simply the action of an external 
electric field on a moving spin. 

If the potential V{t) has spherical symmetric, i.e. 
V(y) = V{r), we have W{r) = ^£y{r). Then the 
general spin-orbit interaction of Eq.|5jl reduces to the 
following familiar form: 



Hso — 



1 1 d 



2m2c^ r dr 



V{r)s»\ 



(6) 



where the orbital angular momentum operator is 1 = r x 
p. In fact, Eq.® is the well-known Thomas SO coupling. 

When our device is made of a two-dimensional electron 
gas (2DEG) in which electrons are strongly confined in 
the y direction by a confining potential V{y)^ such that 
% » %^ then VV{v) « and the electric field 
is almost along the y direction. Furthermore, if V{y) is 
asymmetric with respect to the refiection point y = 0, 
then the matrix element {y)\-§ijV {y)\^ {y)) ^ where 
^'(y) is the basic bound states in the y direction. Un- 
der these conditions, the general SO interaction Eq.Q 
reduces to the Rashba SO interaction formiSliS^ 
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[a{(7 x p) + ((7 X p)a] 



(7) 



where a ^ (^(2/)|^^(2/)|^'(2/)} is the interaction coef- 
ficient. Note that an asymmetrical confining potential 
in the direction perpendicular to the 2DEG (the y di- 
rection) is necessary, otherwise a = and there will 
be no Rashba SO interaction. It is worth to mention 
that the Rashba SO interaction strength a can be tuned 
in an experiment by an external electric field or gate 
voltage, which have already been done in some recent 



A. Without SO interaction 

The second quantized form for the Hamiltonian 
of the Metal-QD-Metal device with non-magnetic leads 
(M = 0), in zero magnetic field (B = 0), and without 
SO interaction (a = 0) can be approximately written in 
the standard Anderson model: 



(8) 



where Hqd is the Hamiltonian for the QD region; Hp is 
for the leads and Ht the coupling between the leads and 
the QD: 

,ns' '^ns'^ms' (9) 

Hp ^^ekpalp^ttkps (10) 



Ht 



H.c 



kf3,n.s 



(11) 



The quantity n„s = diisdns', s =T,i (or -I- and -) is 
the spin index which also describes the spin states, with 
\s) = (1,0)"^ and (0,1)"^ for the spin-up and spin-down 
states, respectively, n is quantum number for eigen- 
states of the single-particle Hamiltonian Hg (Eq.(|3J)) in 
the isolated QD region with eigen-energy e„ — {n\Hs\n); 
k(3 is the quantum index for lead (3 with eigen-energy 
ekp = {k(3\Hs\kf3) {f3 = L,R). tkp,n = {kf}\Hs\n) de- 
scribes coupling strength between the leads and the QD 

region. Quantity Uns,ns' = (ns, 'ms'\ ^^^^_^^yi \ms' , ns) is 
the matrix element for the two-body e-e interaction. Here 
the e-e interaction in the leads region has been neglected. 
Note that when M = B = a = 0, ffs does not depend 
on the Pauli matrix ct, therefore [a, Hs\ = and e„, ekp, 
and ifc/3,n are all independent of the spin-index s. 

For ferromagnetic leads M.p ^ and Hg of Eq.lO con- 
tains a term a • 'M.p. Let's assume has a constant 
value in each lead /3 although ^ M^; in general. 
By calculating matrix elements {kf3s\a • Mplk'^s') = 
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5kk'{^ * the second quantization for this term 

can be easily obtained. Hp becomes: 



k 



o-kpi 



Due to the existence oi &• M^, the state \kf3s) is usuaUy 
not an eigenstate of isolated lead /3, and {aXpgakps) is not 
equal to the Fermi distribution function fp{eki3) even in 
equilibrium. In order to conveniently solve the transport 
problem, we diagonalize Hp above by a rotational unitary 
transformation. This is accomplished by setting: 



a-kfjv \ _ 
a-kpi' 



-sm-j-e^'^fl'^ cos^e-'-^f/^ / \ akpi 



a-kf}]' 



where 9p and <j)p are the directional angles of the FM mo- 
ment M^. Under this rotational transformation, the to- 
tal second quantized Hamiltonian of the Metal-QD-Metal 
device becomes: 

ns,ms' (ns^ms') 

Hf3 = 51 (^'^'^ + sMp)alpgakf)s (12) 



Ht= (cos y 



-ssmyaj^^j I X 



H.c 



where Mp — |M^|. This form of the Hamiltonian has 
been used beforc"^^, but two important comments need 
to be made: (i) In Hamiltonian H12() . the states \k(3s) are 
eigenstates of Hp for isolated leads, hence in equilibrium 

{^ki3s"'ki3s) = fpi^kps) with tkps = efc/3 + sMp. (ii) Af- 
ter the rotational transformation, the spin-up direction 
in the left FM lead, the QD and the right FM lead, are 
all different although they are all aligned in their local 
z-direction. These local coordinate systems are shown in 
Fig. 2. In the QD, the spin-up direction is still in the orig- 
inal z axis; but in the left /right FM leads, the spin-up 
direction {i.e. local z-direction) is aligned with the FM 
moment M^/^ (see Fig. 2). Although this difference of 
spin-up alignment is not important when the QD bridg- 
ing the leads has only a single connection (such as Fig.l), 
it is important if the QD region has double or more con- 
nections (such as Fig. 5). 



B. Rashba SO interaction (I) 

In this and the next subsections, we second quantize 
the Rashba SO interaction which is a major content of 
this paper. The Rashba SO interaction (TJ can be splited 



into two terms: 
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[a{x){(T X p) + ((7 X Yi)a{x)] 

— [a[x)(TzPx + azPxa{x)\ 

Hri + H 



R2 



(13) 



For transport direction along the x axis as shown in Fig.l, 
these two terms have some essential differences. The first 
term Hm gives rise to a spin precession^ while the sec- 
ond term Hji2 does not. In particular, Hm includes a 
(5-function factor at the Metal-QD contacts {x = x^/r, 
see Fig.l) with a = and the QD region with a ^ 02^. 
For this reason it cannot be second quantized by simply 
calculating the matrix element (nsjiJ^ji |ms'). To over- 
come this difficulty one has to choose a new basis set in 
the QD. This will be accomplished in this subsection and 
the Hji2 term will be studied in the next subsection. 

For clarity, the real space single particle Hamiltonian 
considered in this subsection is: 
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V[v) -f CT . M(r) + Hri . (14) 



This is just Eq.Q but with only the Hjn part of the SO 
interactions. We make an unitary transformation with 
the following unitary matrix: 

{1 X < XL 

exp{-iaz J^^ kR(x)dx] xl < x < xr (15) 
exp{-iaz J^" kR{x)dx} xr < x 

where kR{x) = a{x)m* /h^. Here a{x) is permitted to 
have a dependence on spatial coordinate x inside the 
QD, and it is zero outside {xr < x or x < xl, see 
Fig.l). Under this unitary transformation, the original 
basis functions in the QD region, \n f) ~ (^„(r)(l, 0)"^ 
and \n |) = (r)(0, 1)"^, arc transformed to 



\n T)' - u{x)\n T) = e-'-f^. '=«(-)'^>„(r) ^ j (16) 

In i)' = u{x)\n i) = e^'^^\ '=H(-)'i-<^j;(r) « ^ (17) 

These new basis functions are used to second quantize 
((TUl . After the unitary transformation. Hi of Ea. l(n|) 
becomes: 



1 



Hi' = u{x)^Hlu(x) 



Px+Pz 

2m* 



V{r) 



h^k\{x) 

2m* 



cr«M'(r) (18) 



where M'^ = and |M^| = |M_r|, but the directional 
angles of are changed to {Or, 4>r~ 20so) with (j)so = 
J^^ kR{x)dx. 

The essence of the above unitary transformation are 
the following, (i) It is equivalent to choosing a space- 
dependent spin coordinate as shown in Fig. 3a, in which 
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the spin z direction is fixed everywhere but the spin x 
and y directions are dependent on the space position r. 
In different positions along the x axis, the directions of 
the spin x,y axis are rotated. In other words, the uni- 
tary transform changes us to a rotating frame. It is weU- 
known that for an electron moving along the x direction, 
the Rashba term Hri gives rise to a spin precessionii^. 
That is, the spin component in the x-y plane will rotate 
as the electron moves along the x direction, therefore 
the electron spin is usually not invariant. However, in 
the rotating frame which follows the spin precession, the 
spin is invariant hence [H^ , ax/y/z] = is satisfied in the 
QD region, (ii) The Rashba interaction Hm can cause 
an energy split between spin-up and spin-down states for 
non-zero fc^, as shown by the energy dispersion in the 
left panel of Fig.SbSS*^. The above unitary transforma- 
tion recovers the alignment of the two dispersion curves 
so that the right panel of Fig. 3b is obtained. Therefore, 
after the unitary transformation, the new Hamilton 
appears to be completely the same as the Hamiltonian 
without the Rashba interaction except a rotation 

of the magnetic moment M/j and a potential energy dif- 
ference — ^ 2m*^'' ~ — m*a^(x)/(2fi,^) which is a simple 
constant if a{x) is independent of x. Using the same 
method as that of the last subsection, the second quan- 
tization of Ea. H18() is easily obtained: 

H — Hqo + Hp + Ht 

I3=L,R. 

n,s 

= "^i^kf! + sM0)alp^aki3s (19) 

k.s 



k,7l,S 



tkLnicos y'^Ils - 5 sin %aL,-)e'"'^^/2(i„ 

t ■ t \ 

2 -kRs - ssm— a[^j) 



+ t.H„(cos ^4^, - s sin ^ai^,)e-^-/'e-*-d„, 

+H.c] (20) 



This is one of the main results of this paper. Here the 
Rashba interaction Hj^i gives rise to an extra spin depen- 
dent phase factor —s4>so in l|2()|l : it is —4>so for s =], and 
is +(t>so for s —I. Note that the term with this phase fac- 
tor satisfies the time-reversal invariance while = 0, 
i.e. [T,Ht\ = where T is the time- reversal operator 
This is an expected property because the Rashba SO in- 
teraction in real space (Eq.©) does satisfy time- reversal 
invariance (see the Appendix). We emphasize that the 
phase factor —s(j)so in Ea. (|20|l is fundamentally different 
from the phase factor caused by magnetic flux in systems 
such as the AB-ring: the latter is independent of spin s 
and it destroys time-reversal symmetry. 

For the special case where fcfl;(a;)=fci{=constant, i.e. 
independent of coordinate x of the scattering region, 
we have (f)so = kj^ x (a;^ — x^). Then, re-defining 
-iskRXL^ _^ Pgg^j space Hamiltonian (|20|l can 



be rewritten in a symmetric manner: 

n,s 



(21) 



ifc/3n(c0S Y"-kl3s - S Sm ycj^^j) 



Ht = 



C. Rashba SO interaction (II) 



Now we second quantize the second term of the 
Rashba interaction Ea. (|13|l . Hfj2 = —SdElEsS^^ which 
can be accomplished by calculating the matrix elements 
{ms'\u{xy Hji2u{x)\ns) = {ms'\H'j^2\^s)- If = s, this 
matrix element is exactly zero. Hence we only need to 
calculate the non-diagonal matrix elements, and they are: 



{ml\H'j,2\nr) 
m* 

{nl\H'j,2\m^) 
—hkfi 



dre-^^^-->™(r)p,<^„(r)^Cn (22) 



TO' 
TO* 

{n UH'n2\m i) ^ t 
{mUH'R2\ni) = -tZ 



dre 2*'=fi^(p„(r)p^^„(r) 
dre-2''=«>„(r)p,^„(r) = (23) 



±SO* 

7nn 



(24) 
(25) 



Here (as well as in below) we have assumed that kR{x) 
(or a{x)) to be independent of x, but even if kji{x) de- 
pends on x, all results are completely the same. With 
the above matrix elements (I22l25f) . the second quantized 
form of Hii2 is: 



Hr2 - ^ i^, 

rn.n{ra<iri) 



4a '^"T 



H.c. 



which can be written in a more compact form: 

iffl2=^CXi^«T+^^-C- (26) 



where it is important to realize that = — 

Ea. (|25|) is another main result of this paper. 

Some general characteristics of Ea. H26() are in order, (i) 
The property = — i^'^ for the matrix elements origi- 
nate from the time-reversal invariance of the original real 
space Rashba Hamiltoniani^ Using this property, we can 
exactly prove that the second quantized form of Hr2, 
Ea. (|26|) . indeed satisfies the time-reversal invariancei^ 
(ii) If n = m, we have = — t^'^ hence i^"^ must van- 
ish. This means that the Rashba SO interaction can not 
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induce any intra-level spin flip, i.e. it can not give rise to 
a transition (n |) — > (n J,) in which the level index n is the 
same. Therefore, the SO interaction (|26|l is fundamen- 
tally different from that of an external magnetic field: 
magnetic field can cause intra-level spin-flip, provide a 
Zeeman energy that relieves spin-degeneracy, and induce 
spin-polarization in an isolated QD (see subsection II-E 
below), (iii) The Rashba SO interaction H26|l can cause 
spin flips between different energy levels. This inter-level 
spin-flip coupling is similar to the intersubband mixing 
in real space which has been studied in previous worksi^. 
Despite of the inter-level spin flips, the system is still 
at least two-fold degenerate for any eigenstates because 
^mn = ~^nm- This guarantees that at equilibrium an 
isolated QD has no spin polarization. In the Appendix, 
the general properties of Ea. (|26|) are further discussed, 
(iv) In fact, because all spin-orbit couplings satisfy time- 
reversal invariance, the above properties and matrix ele- 
ments tmn = —tnm must hold truc in general. In this re- 
gard, we note that there existed papers where tmn = tnm 
and intra-level spin-flips are allowed: these effects cannot 
come from SO interactions as some times claimed. 

Let's estimate the value of t^^. Consider a square 
QD with linear size W. The eigenstates are ipn{^) = 
■j|r sin sin , hence can be easily calculated 
from Eas. 1221) . For parameters W = lOOnm, a = 3 x 
lO^^^eVm and m* = 0.036me, the intradot level spacing 



A. 



2m* 

rough estimate of |<f°j 



ImeV. This is to be compared with a 



D. Electron-electron Coulomb interaction 

In order to second quantize the Rashba SO interaction, 
we have introduced an unitary transformation defined by 
Eas. H15|l . Does this transformation affect the familiar 
second quantized form of the e-e interaction? Here we 
show it does not. 

Starting from the two-body e-e interaction in real 
space: 



Hi{vi,V2,...Vn)^ X! 



2\vi -r,|2 



we apply the unitary transformation and the new Hamil- 
tonian H'j is: 



- y 



u{xj)u{xi) 



interaction in its familiar form, 



(27) 



where the matrix element Uns,ms' is: 

^2 



Un 



{ns, ms 



2 ri -ra 



\ms', ns) 



E. External magnetic field 

The unitary transformations Eas. (|15|l does affect the 
second quantized form of the external magnetic fields 
(T • B. Consider an arbitrary external magnetic field 
B = {B^, By, Bz) where B^jyj^ is projection in the xjyjz 
direction. 

First, we investigate the z-direction element B^. Un- 
der the unitary transformation Eas. H15|l . the term a^Bz 
changes to: 



which means dzBz does not change under the unitary 
transformation. Therefore its second quantized form is 
still 



YsBzdlj„ 



(28) 



Second, we investigate the x-direction element B^, i.e 
the term cr^B^ in the Hamiltonian. After the unitary 
transformation Eas. (|15|l . {x){(JxBx)u{x) ^ ^xBx so 
that it is affected by the transformation. The matrix 
elements {nis' \u^ {x)axBxu{x)\ns) are found to be: 



(to t \u^axBxu\n t) 
(to I \u^axBxu\n t) 

(n I \u^(7xBxu\m t) 
{n t \u^(7xBxu\m |) 



(to i \u^axBxu\n j) = (29) 
dre-'^''''^^„,{vMr)Bx 

B 
mr 
B 



4-a R 



4-a n 

(to t \u'^crxBxu\n [) 

iB* r> 

''mn^x 



(30) 
(31) 



(32) 



Hence, the second quantized form of cTxBx is: 



inn 



J2 Bxt 

m,n(n<m) 

+ J2Bxt^ndUdn^+H.c. 

n 

or it can be written in a more compact form, 
J^BxCndl^dn^+H.c. 



(33) 



This means H'j = Hj and the unitary transformation 
does not affect the form of the e-e interaction. We there- 
fore can directly write down the second quantized e-e 



with = if„j. Note that the form of Eq.(|2SIl is very 
similar to the Rashba term Hii2, Ea. H26|) . but there ex- 
ists an essential difference. Namely for magnetic field 



in the x-direction, = in Ea. H33|) : while for the 
Rashba term, = — ifj"^ in Ea. l(^ . We emphasize 
that this is an essential difference because of two reasons: 
(i) The magnetic field term destroys time-reversal invari- 
ance, it provides a Zeeman energy that breaks the spin 
degeneracy of the energy levels, and it can induce a spin- 
polarization in equilibrium. In contrary, the Rashba term 
Hii2 satisfies time-reversal invariance and maintains the 
two degeneracies, (ii) When n = m, can be non-zero 
so that intra-level spin flips are possible. Furthermore, 
the term is usually the largest term in the sum of 
Ea. (|S^ . e.g. = Smn at = 0. But for Rashba 

interaction (|26|l . tf^„ must vanish as discussed before so 
that it cannot cause intra-level spin flip. We therefore 
comment that interactions of the following form, which 
has been used in some previous literature. 



where 



ns^ms' {ns^ms') 



+ 



(36) 



k,s 



+ / t • t ^ 

tkf3n{cOS yKps - S Sm y a^^j) 



(37) 



tdld 



does not represent SO interaction. Rather, it describes a 
magnetic field pointing to the x-direction. 

In order to estimate the value of t^,^, we consider a 
rectangular QD with length L and width W. can be 
obtained as: 



dxe 



~2ikRLx 



sm nixTrx sin UxTtx 



where — —t^"^ and = This Hamiltonian 

is the central result of this paper. The Rashba SO in- 
teraction causes two effects, (i) It gives rise to an extra 
phase factor —skjixp in the hopping matrix element be- 
tween the leads and the QD. Note that this phase factor is 
dependent on the electronic spin s, and it is essential dif- 
ferent from the usual phase factor due to a magnetic flux 
which is independent of s. (ii) The Rashba SO interaction 
causes an inter-level spin-flip term with the strength t*°„ 
and it cannot cause intra-level spin flips. Time- reversal 
invariance is maintained by the SO interaction which is 
essential different from the effect of an external magnetic 
field. 



Fig. 4 plots numerical results for versus parameter 
fcfli obtained this way. As UrL increases, more are in 
action. If parameters a = 2*lQ~^^eVm and L = lOOnm, 
we have kjf L « 1 for m* ~ 0.036me. For this kuL value, 
only tfi and are significant. 

Finally, the second quantization of the By term is com- 
pletely the same as that for the x-direction, hence its 
second quantized form is the same as Eg . 1)3311 . 



F. Brief summary 

Collecting all the pieces of second quantization which 
we have carried out in this section, for a device in the 
form of Metal-QD-Metal where there exists Rashba SO 
and e-e interactions in the QD, the Metal leads are mag- 
netic material, and there exists an external magnetic field 
B, Hamiltonian becomes: 

H = Hqd+ E Hp + Ht (34) 



IV. EXAMPLE: TRANSPORT PROPERTIES OF 
AN AB RING WITH RASHBA SO 
INTERACTION 

As an example of applying the second quantized Hamil- 
tonian Eas. H34l37|) . we now investigate quantum trans- 
port properties of a modified AB ring shown in Fig. 5. A 
QD sits on one arm of the ring and Rashba SO interac- 
tion exists inside the QD. No Rashba interaction exists 
on the other arm of the ring. The ring is connected to 
the outside world by two normal metal leads. AB rings 
with an embedded QD have been studied by many pre- 
vious works^Si^LSSi^. Some interesting phenomena such 
as Fano resonance22i^2i2S have been discovered in such a 
device. The effect of Rashba interaction has not been 
studied so far, and we have found that it leads to in- 
teresting transport behavior. In particular, a substan- 
tial spin-polarized current or conductance is induced by 
a combined effect of the Rashba SO interaction and a 
magnetic flux cj) threading through the AB ring. The 
direction of spin-polarization and its strength are easily 
controllable by or by a gate voltage. 

The Hamiltonian of our AB ring (Fig. 5) can be written 
using various pieces of the general Hamiltonian Eas. (|34t - 
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H 



k,s,0{l3=L,R) s 



k,s 



H.c.{38) 



As discussed previously, the physical meaning of each 
term is clear: the first term describes the normal metal 
leads; the second term is for the QD which has a single 
energy level with spin index s; the third term is the in- 
tradot e-e Coulomb interaction with a constant strength 
[/; the fourth term is for the arm of the ring without 
the QD; and the fifth term is the coupling between the 
leads and the QD. Due to the Rashba SO interaction, ac- 
cording to Eq. (|37() there is a spin-dependent phase factor 
—skuL in the hopping matrix element tnd on the fifth 
term. Since we only consider one level in the QD, the 
inter- level spin- flip term of Ea. (|35|l does not appear here. 
This is equivalent to neglecting the intersubband mixing 
as in some previous work*". We emphasize that both the 
e-e Coulomb interaction and Rashba SO interaction are 
considered, different from previous studies of the Rashba 
SO interaction where e-e interaction was neglected. In- 
deed, our second quantized Hamiltonian in the spectral 
space Eas. (|34l37f) allows us to consider both effects to- 
gether. Finally, the magnetic flux $ threading the AB 
ring gives rise to a familiar spin-independent phase factor 
(j) = 27r4>/$o in the matrix element tud- 

The quantum transport problem described by Hamil- 
tonian H38|) can be solved by standard many-body tech- 
niques. In the following we calculate charge current using 
the standard Keldysh non-equilibrium Green's function 
method. Following Ref l23. the charge current flowing 
from the left lead into the AB ring, contributed by spin- 
up or spin-down electrons, can be derived as: 

= ¥ / [^^^^^dLsi^'^ + ti^RG<M] (39) 

where Keldysh Green's function G'^{uj) is the Fourier 
transform of G^{t) whose definition is: 

k' k 

G^dsit) ^ ^{dliO)Y,akps{t)) 

k 

G^dAt) ^ ^dlmsit)) . (40) 

To solve , we first calculate the retarded Green func- 
tions using the Dyson equation: 

G: = g:+g:s:G: (4i) 

and the Green's function Gl is a 3 x 3 matrix defined as: 



G 

g: = I G 



LLs 
r 

RLs 
^dLs 



^LRs 
^RRs 
^dRs 



Lds 
'-'Rds 
^dds 



(42) 



In Ea. (|41|l . is the Green's function of the system with- 
out coupling between the leads and the QD (i.e. when 
tj^R = thd = tRd = 0). It can be obtained exactly as: 



-inp 
:{uj) = I -iTrp 

g-,,M 



(43) 



where q^jjAlo) = " '^\/^'^^"'',t\ and ric is the intradot 
electron occupation number at state s. p in Ea. H43|l is 
the density of state of the leads. The self-energy Sj(ijj) 
in Eq.gH) isi^ 



tLR thd 
tRds 



f* n 
f* f* n 

''Ld ''Rds " 



(44) 



where inds = tude-^'''^^ e^'l' . Using Eas. iPEll . G^ can 
easy be obtained by solving the Dyson's equation H41() 
as: G^ = (gr^-S-)-i. 

After solving Gs(w), the Keldysh Green's function 
G^(ijj) can be obtained straightforwardly from the stan- 
dard Keldysh equation: 



G< = (1 + G;s;)g<(l + S:G:) + G:S<G: 

= G%-^g<griG: + G:S<G: (45) 



For our present case, — and g^ g^ Ss diag- 
onal with ^'p'pl^pps^'Bl = 2i/^(tj)/7rp {(i = L,R) and 
^Mh^ds^dll = 0> where fp{u^) = l/(e(-M.)/'=i.T + 
is the Fermi distribution function in lead j3. As the last 
step, the intradot electron occupation number needs to 
be solved self-consistently with the self-consistent equa- 
tionn, = -z/|^G<,(^). 

In the following we present numerical results. Fig. 6 
shows the total linear conductance G = ^ ver- 

dV 

sus the intradot level position at zero magnetic flux 
{(f) = 0) but with different Rashba interaction strength 
kuL: fcflL = (solid), 7r/4 (dashed), 7r/2 (dotted), 
37r/4 (dash-dotted), and tt (dash-dot-dotted). The curves 
are dominated by two Coulomb peaks at = and 
Erf = —U. When there is no Rashba SO interaction {i.e. 
kuL = 0, solid curve), these two peaks show a typical 
Fano resonance shape due to the interference of electrons 
passing the two arms of the AB ring, in agreement with 
previous theoretical and experimental studies^SiSi. It is 
interesting to discover that this Fano resonance can be 
strongly affected by the Rashba SO interaction. Increas- 
ing the Rashba parameter kpiL from 0, the Fano reso- 
nance is decreased and it can completely disappear at 
kfiL = 7r/2 (dotted curve). Further increasing fc^jL, the 
Fano resonance rises up again but with an oppositive 
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Fano factor, for example at kj^L = tt (dash-dot-dotted 
curve). 

In order to understand these results, we investigate the 
interference term of total transmission probability, which 
is approximatively proportional to ^ J^s cos{A0 +skiiL). 
Here A9 is the phase difference of the transmission am- 
plitude through the two arms, and it varies from to 
tt/2 and finally to tt as is moved from —oo to and 
finally to oo. This clearly shows that the total trans- 
mission probability is indeed having a Fano asymmetric 
resonance shape when kjiL = or tt. On the other hand, 
it is symmetric at k^L = tt/2 or 3tt/2. Hence the Rashba 
SO interaction can alter the Fano resonance shape in sub- 
stantial ways. 

Next, in the three panels of Fig. 7 and Fig. 8, we 
plot conductance Gg and spin polarization rj = {G^ — 
Gi)/ (G| -I- Gi) versus magnetic flux for three values of 
ed — 1, 0, —1, respectively. These values of are near the 
right Coulomb peak of Fig. 6. In Fig. 7, the thick curves 
are for G| and thin curves for Gi, the solid, dashed, and 
dotted curves correspond different values of the Rashba 
parameter kuL = 0, 7r/4, and tt/2. Figs. 7, 8 clearly show 
that if cither one of the two parameters {(j) and k^/L) 
vanishes, the transport current has no spin polarization 
so that 77 = and G| = G|. However, when both pa- 
rameters are nonzero, a substantial spin-polarized con- 
ductance is found and rj can be as large as 90% for the 
given set of system parameters (Figs. 8a, c). 

Importantly, in the present device the direction of spin- 
polarization and its strength (77) are easily controllable by 
varying system parameters which are experimentally ac- 
cessible, (i) By varying the magnetic flux 0: when is 
tuned from —tt/2 to 7r/2 (or from tt/2 to 3tt/2), the polar- 
ization 77 strongly varies from a large positive value to a 
negative value or vice versa, (ii) By varying the intradot 
level €(1 using a gate voltage: when is moved from one 
side to another side of a Coulomb peak, the polarization 
T] can be tuned from its largest positive value to its largest 
negative value or vice versa. Numerically we found that 
one only needs to change ed by a small amount to see the 
polarization change, namely a few half-widths T of the 
Coulomb peak (parameters used in Figs. 7, 8 correspond 
to F = 2TTp\tpd\'^ ~ !)• This means that in an experi- 
ment one only needs to slightly vary the gate voltage to 
change ry from 1 to —1 or vice versa. Furthermore, we 
note that when polarization rj reaches its largest value, 
the conductance itself is still large, e.g. Gf or G| can 
exceed over 0.8e^//i (see Fig. 7a). 

Finally, we estimate if the parameter k^L can reach a 
value of tt/2 in the present experimental technology so 
that the above theoretical predictions can be observed 
experimentally. Assuming the Rashba SO interaction 
strength a ~ 3 x lO^^^eVm, which is the reported 
value for some semiconductorsiSiSSiffi, kjf ~ m*a/h^ k, 
0.015/nm for m* = 0.036me. Then, if the length of the 
QD is the typical value lOOnm, k^L k, 1.5. Therefore 
we conclude that k^L can reach a value ^ tt/2 or larger 
experimentally. 



V. CONCLUSION 

In this paper we have derived a second quantized 
Hamiltonian in spectral space for a general device struc- 
ture of Metal-QD-Metal configuration, including the 
spin-orbital and e-e interactions. In other words, we ex- 
tended the standard Anderson Hamiltonian to the case 
where the central device region has a (Rashba) spin- 
orbital interaction. We discovered that the Rashba SO 
interaction causes two changes: (i) It gives rise to an ex- 
tra spin dependent phase factor —skuxp in the coupling 
matrix elements between the leads and the quantum dot. 
(ii) The Rashba SO interaction causes an inter-level spin- 
flip term with strength = — i^^, and it cannot cause 
any intra- level spin flips. 

The spectral form of the Hamiltonian is very important 
as it allows the analysis of many complicated quantum 
transport problems involving SO and e-e interactions, by 
using the well established many-body Green's function 
theoretical techniques such as the Keldysh nonequilib- 
rium Green's function formalism. On the other hand, it 
would be much more difficult to carry out similar inves- 
tigations using a real space Hamiltonian, especially if e-e 
interactions are present. 

As an example, we investigated quantum transport 
properties of a AB ring in which a QD having Rashba 
SO and e-e interactions is embedded in one arm of the 
ring. A substantial spin-polarized current or conductance 
emerges in this device due to a combined effect of a mag- 
netic flux and the Rashba SO interaction. In particular, 
the direction of the spin polarization and the strength rj 
can be easily controlled by a number of experimentally 
accessible parameters. 
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Appendix 

In this appendix, we collect some general properties 
of the spin-orbit coupling. Although these properties 
should be well-known beforc"'^^-, we believe it is useful to 
put them in a form that is easily accessible. In addi- 
tion, these properties hold for all spin-orbit interactions 
including the Rashba SO interaction. 

(1). The SO interaction Hamiltonian Hgo^ Eq.(|SJl, 
satisfies time-reversal invariance. In other words, Hgo 
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commutes with the time-reversal operator T = —iayK 
(where K is the complex-conjugation operator). Indeed, 
using the Hamiltonian Hgo of Eq. © , it is easy to prove 

(2). When a system has spin-orbit coupling, each 
eigenenergy level is still at least two- fold degenerate, i.e. 
the so-called Kramers degeneracy exists. Briefly this can 
be proved as follows. We start from the Hamiltonian H 
of Eqs.p but setting B = M(r) = in assume 
(^ji(ri,r2, ...rjv) is an eigen-state of H so that H\ipn) — 
En\^n)- Since H is time- reversal invariant [TH = HT), 
we have HT\ipn) — TH\ipn) — EnT\(pn). Hence state 
T\ipn) is also an eigen-state with the same eigen en- 
ergy En as that of the state \ipn). Furthermore, one has 

{ipn\T\ipn) = {ipn,Tipn) = {T'^ipn,Tip„) = - {ip„\T\ipri) , 

hence {Lpn\T\ipn) = 0. This means state T\ipn) is orthog- 
onal to \ipn). Therefore, although spin is no longer a good 
quantum number when SO interaction exists, the system 
is still at least two-fold degenerate for any of its eigen 
state. 



(3). At equilibrium, any spin-orbit coupling cannot 
induce a spontaneous spin polarization. We prove this as 
follows. Since the system is in equilibrium, the two- fold 
degenerate eigen-states \ipn) and T\ipn) have the same 
occupation probability p(_E„). Then, the average of spin 
polarization at an arbitrary direction n can be calculated 
as: 

n 

n 

= '^p{En)[{(Pn\crn\^n) + {T(pn\ - T(7,n |(p„)] 
n 

= ^P{En)[{'fin\crn\'Pn) - iVnlcTnlVn)] 
n 

= . 

Therefore, at equilibrium no spin-orbit coupling can in- 
duce a spontaneous spin polarization in any direction. 
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This is a good approximation for temperatures higher than 

the Kondo temperature. 

FIG. 1: (a) Schematic diagram for a Metal-QD-Metal device 
configuration where the QD is weakly coupled to two fer- 
romagnetic leads, (b) Schematic diagram for the scattering 
potential along the x-direction. The Rashba SO interaction 

is assumed to only exist in the central QD region, i.e. a = 
in for regions with x < xl and x > xr. 



FIG. 2: Schematic diagram for the spin coordinates, i.e. the 
spin-up direction in the left lead, the center region, and the 
right lead, respectively. 



FIG. 3: (a). Schematic diagram for the spin coordinate axis 
in different position. Hero the x and y spin directions arc ro- 
tated along the a:-axis in space, (b). Schematic diagram for 
the dispersion relation before and after the unitary transfor- 
mation. 
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FIG. 4: The spin-flip coupling strength versus the Rashba 
SO interaction UrL. As knL increases, more modes (m) are 
playing a role. 



FIG. 5: Schematic diagram for the modified AB ring device: 
two normal leads arc couplied to the center ring, a magnetic 
flux threads the the ring, and a QD is embedded to one arm 
of the ring. 



FIG. 6: Conductance G versus the intra-dot level td for 
UrL = (the black solid curve), 7r/4 (the red dashed curve), 
7r/2 (the magenta dotted curve), 37r/4 (the blue dash-dotted 
curve), and tt (the purple dash-dot-dotted curve), respec- 
tively. Other parameters are: tad = tLd = 0.4, Ilr = 0.1, 
pi^ = Pr = 1, ksT = 0.0001, U = 5, and 4> = 0. 



FIG. 7: Conductance G| and G| versus magnetic flux (j) for 
several intradot levels: (a) = 1; (b) td = 0; (c) = — 1. 
The solid, dashed, and dotted curves correspond to feijL = 0, 
7r/4, and 7r/2, respectively. The thick curves are and thin 
curves are Gj. Other parameters are the same as those of 
Fig.6. 



FIG. 8: Spin polarization 77 versus magnetic flux </> for fc^L — 
(dashed curve), n/4 (solid curves), and 7r/2 (dotted curve), 
respectively. Other parameters are the same as those of Fig.6. 
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